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I. INTRODUCTION 



In this paper we use an exactly solvable model of disorder to compare different theoret- 
ical approaches used to investigate the critical behaviour of disordered (possibly electronic) 
systems. A prominent example of a disordered critical point, and a source of many challeng- 
ing problems, is that governing the plateau transitions in the integer quantum Hall effect 
(IQHE) — see for example J}]]. Despite the numerous field theoretic approaches to this 
problem based on either Pruisken's non-linear sigma model [0-0] or the Chalker-Coddington 
network model || and its superspin chain descendants |6|-[I0| , the nature of the critical point 
remains elusive and inaccessible by perturbation theory [|TT],|T2[. This is to be contrasted 
with the Ising model with weak bond disorder [|l3Hi5|], for example, where the disorder 
renormalizes to zero and may be treated perturbatively. In order to gain insight into such 
problems it is imperative to fully exploit, and indeed develop, the methods available for 
treating disordered systems, and to obtain as many non-perturbative results as possible. 
In low-dimensional systems we may be aided in this pursuit by the availability of powerful 
techniques such as conformal field theory [|i6| --|i"8|l and the Bethe ansatz fl9l,[20|. 

Recent studies of disordered critical points have included Dirac fermions subjected to 
various random potentials [f21~H 24|l , the random XY model [25|, the plateau transitions oc- 
curring in the spin quantum Hall effect and its relation to critical percolation P^-EHf, and 



the Nishimori line in the random bond Ising model [29]. In all of these examples the critical 



points are of a non-perturbative nature and are described by non-trivial field theories. Our 
knowledge of the critical behaviour in these theories differs somewhat, and the further elu- 
cidation of their detailed properties is a valuable enterprise. Indeed, in problems of disorder 
related to the localization of quantum particles, one needs to depart from the critical point 
in order to calculate the diffusion propagator or (for problems with a singular density of 
states) the energy dependence of the density of states; at present this program has been 



successfully implemented only for the random XY model |25 



We turn our attention now to the problem at the heart of this paper, namely two- 
dimensional (Euclidean) Dirac fermions in a random non-Abelian gauge potential. This 
model has the virtue of being amenable to a variety of non-perturbative approaches and was 
originally introduced in Abelian form as part of an attempt to describe the plateau transi- 
tions in the IQHE |^TJ. The non-Abelian version of the problem appeared in a treatment of 
disordered d-wave superconductivity [|30| -|32"H and has been the subject of deeper investigation 
and refinement [[33H 4011. Building on the foundations of these previous studies, we address 
this problem by means of three independent non-perturbative field-theoretic techniques: 
these are based upon the commonly used replica |^T| and supersymmetric methods [[42||43| 
together with a third (model specific) approach valid in the limit of strong disorder [ 39"| . In 
our study of the four-point correlation functions of the local density of states we are able 
to demonstrate that the three methods yield coincident results. We note that whilst the 
reliability of the replica approach outside of perturbation theory has been frequently and le- 
gitimately questioned — see for example [^-[£5| — we use it here (in a quite straightforward 
manner) to reproduce the results of the (mathematically more rigorous) supersymmetric 
approach. Although the implementation of the replica method remains a delicate issue in 
general, and in many cases one must adopt various (replica) symmetry breaking schemes - 
see references [[K|^7j for recent examples in random matrix theory — it is noteworthy that 
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the disordered non-Abelian Dirac fermion problem renders itself here to a textbook (replica) 
treatment — albeit with the benefit of hindsight. 

The structure of this paper is as follows: in section we provide a brief outline of 
the disordered non-Abelian Dirac fermion problem. In section |T| we focus our attention 
on the four-point correlation functions of the local density of states (LDOS) as encoded 
in the so-called Q-matrix. This section is subdivided into three main subsections, each of 
which is devoted to a different non-perturbative field theoretical approach to the disordered 
Dirac fermion problem; subsection [III A| deals with the replica approach, [IIIB| deals with 
the so-called strong disorder approach and [111 (J| deals with the supersymmetric approach. 
Within each of these subsections we provide an outline of the theoretical approach and 
the form of the appropriate Q-matrix, together with a discussion of the relevant conformal 
dimensions and four-point correlation functions. In order to increase the transparency of 
these subsections we have relegated many of the important (but arguably involved) technical 
details on the solution of the relevant WZNW models into a rather substantial appendix. 
The interested reader will find gathered in this appendix the solutions to the Knizhnik- 
Zamolodchikov equations which arise in this work. In section [IV] we continue our study of 
the supersymmetric approach and provide a remarkably simple free field representation of 
the disorder averaged theory comprising of a two-component symplectic fermion and a pair 
of free bosons (one compact and the other non-compact). The two-component symplectic 
fermion with central charge c = —2 arises in the theoretical description of dense polymers and 
we demonstrate that the twist operators of the latter theory play a fundamental role in the 
non-Abelian Dirac fermion problem. In section |V| we discuss the convergence of approaches 
to the disordered non-Abelian Dirac fermion problem. Finally we present concluding remarks 
and technical appendices. In particular, appendix |A| is devoted to the su(N) k WZNW model 
and is subdivided into our considerations of the u(0)k WZNW model — arising in the replica 
approach — and the ~su(N)_2N WZNW model — arising in the strong disorder approach. 
Appendix |B| is devoted to the 6sp(2\2)k WZNW model which arises in the supersymmetric 
approach. 



II. NON-ABELIAN DIRAC FERMIONS 

We consider Nc colours of two-dimensional massless Dirac fermions minimally coupled 
to a non-Abelian gauge field e su(Nc) : 

, N C 

s= d 2 zj2 v^rV 3 , (2.1) 

Q,/3 = l 

where the gauge covariant derivative takes the form 

paf* = ^(5 a(3 d fl -iAf). (2.2) 

The 7 M matrices form a two-dimensional representation of the Clifford algebra {7 M ,7 i/ } = 
2g fiu with Euclidean metric g^ v = diag(l,l), and the gauge field A°f = A^t® 13 may be 
expanded in terms of the generators r a of su(N c ). Physical quantities are obtained by 
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disorder averaging products of Green's functions. We use the distribution functional 

P[A,] oc e~ s ^\ S[A,] = - [ d^TrA,(OK(0 (2-3) 

9 A J 

representing the usual choice of ^-correlated Gaussian white noise for the random vector 
potential. In section [III] we shall compare three different approaches for evaluating the 
disorder averaged correlation functions. 



III. DISORDER AVERAGED CORRELATION FUNCTIONS 



The local density of states (LDOS) of the Dirac fermion problem is given by [3C--32 



N C 

p{v) = Y,[RiL a + LlR a ]. (3.1) 

a=l 

In sections [111 A^PTTC| we shall investigate the disorder averaged multipoint correlation func- 
tions of the LDOS. In particular, we shall focus our attention on the long distance properties 
(obtained for example by integrating out the massive modes) encoded in the so-called Q- 
fields: 

N c N c 

Q~J^RlL a , Qt-J^LlRa (3.2) 

a=l a=l 

We shall study the critical correlation functions of the Q-fields by three different approaches; 
in each case our Q-fields will be represented by Q-matrices governed by appropriate WZNW 
models. 



A. Replica Approach 

In order to average over disorder we introduce Np flavours (or replicas) of the massless 
Dirac fermions appearing in equation (|2.1|) : 

N F N c 

5 Replica = d 2 £j2Yl i^W^iP* (3-3) 
J i=l a, /3=1 

The (replicated) action ( |3.3| ) enjoys a global SU (Np) xSU (Nc) xU(l) symmetry and may be 
recast using Witten's non-Abelian bosonization rules |48|] — for reviews of the bosonization 
approach see for example It is well known that the free Dirac action (in the absence 

any gauge fields) with the flavour-colour symmetry described above may be represented as 
the sum of three critical Wess-Zumino-Novikov-Witten (WZNW) models of the form 

W k [g] = d 2 ZTr'(d»g-%g) + kT[g] (3.4) 
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where (the WZNW term) T reads 

T[9] = 2^J ^^"^'i^gg^dugg-^g). (3.5) 

The celebrated equivalence between the free (flavoured-coloured) Dirac fermion theory and 
the bosonic WZNW models may be written in the symbolic form 

Free Dirac Fermions = su(N c )n f x su(N f ) Nc x u(1) (3.6) 

where it is understood that the chiral blocks of the free Dirac theory are obtained as products 
of the chiral blocks of the WZNW models; here % denotes the group manifold g and the so- 
called level k of the WZNW action ( |3.4| ). That this is indeed an equality between chiral blocks 



is seen most clearly in the work of Fuchs [51|. Of paramount importance in the application 



of ( |3.6| ) to our disordered problem is that the disorder — the su(iVc) gauge potential - 
couples only to (currents from) the su(Nc)n f sector. Averaging over disorder (equivalently 
integrating over the su(Nc) gauge potential) generates a (quadratic) interaction between 
su(iVc) currents only. This interaction scales to the strong coupling regime where a mass 
gap M ~ exp[— 2n/NcX] is dynamically generated in the su(Nc)n f sector of the theory. As 
follows from the decomposition ( |3.6j ) the massless degrees of freedom are SU(iVc) singlets 
and are described by the remaining WZNW models su(N F ) Nc x u(l) ~ u(N F ) Nc . Upon 
integrating out the massive (colour) degrees of freedom, the fermion bilinears are expressed 
in terms of the so-called Q-matrix 

1 N C 

^~lE<A, ? (3-7) 

a=l 

whose indices reside in the flavour (replica) space; this is simply a replicated version of the 
Q-field introduced to describe the LDOS in equations ( |3.1| ) and (|3.2|). The Q-matrix assumes 
values in the group XJ(Np) and is governed by a (critical) effective action of the WZNW 
form (|3.4|) with level k = Nc- The WZNW model thus plays an analogous role to the sigma 



model in the conventional theory of localization [42 - it is an effective action for the slow 



degrees of freedom, once the fast degrees of freedom have been integrated out. 



1. Conformal Dimensions 

In the replica approach, the conformal dimension of the Q-matrix is that of a primary field 
of the (N F — > 0) u(Np)n c WZNW model transforming in the fundamental representation 
— see equation (|A15| ) of appendix [A} 

h Q = ^ (3.8) 



As we shall see in subsections IT I IB 1| and [111 C 1| , the conformal dimension of the Q-matrix 



(LDOS) is given by (|3.8|) in all three treatments of the disorder. This ensures the equality 
of their two-point and three-point correlation functions. It is well known in CFT however, 
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that the four-point correlation functions are not determined solely by scaling dimensions, 
but in general have a non-trivial dependence on the so-called anharmonic ratios (A3) |52 



in 



CFTs of the WZNW form, this dependence is obtained by solving the appropriate Knizhnik- 
Zamolodchikov equations J53|. In order to compare different theoretical approaches to the 
random Dirac fermion problem it is thus essential to study the four-point correlation func- 
tions of the Q-matrix. 



2. Correlation Functions of the Q -field 

The correlation functions of the Q-matrix are those pertaining to the u(Nf)n c WZNW 
model in the replica (N F — > 0) limit, and are thus obtained by solving the u(N F ) Nc 
Knizhnik-Zamolodchikov equations ||53|| . This fact was recognised in the early replica treat- 



ment of the random Dirac fermion problem p0| , |31|| . The approach was subsequently refined 



to accommodate the logarithmic solutions to the Knizhnik-Zamolodchikov equations which 
appear in the replica limit ||33|| . Despite the advances made in these works they suffer (in 



places) from a naive implementation of the replica trick; an extraneous trace over replica 
indices was performed in a number of instances. In order to address these issues, and to 
demonstrate the convergence of approaches to the random Dirac fermion problem, we red- 
erive a number of results and adopt a uniform notation throughout. 

The four-point correlation function of the Q-matrix admits the U(N F ) x U(N F ) invariant 
decomposition |33|,|53| (see appendix [A 1| ): 



(Q rifl (l)Q|; 2f2 (2)Qj3 f3 (3)Q r4f4 (4)) = \z u zk\- 2/n ° £ r^F^z, z) (3.9) 



where r denotes the ordered sequence of flavour (replica) indices ri,r 2 ,r 3 ,r^, and where the 
invariant tensors I\ and I 2 are defined as 7[ = S ri,r2 5 r ' 3 ' ri and 1 2 = S ri,r ' i 5 r " 2,r4 , together with 
similar equations for I\ and I<i- The anharmonic ratio z is defined as z = ^12^34/^14^32 and 
similarly for z. The functions F%j (z, z) are single- valued combinations of the solutions to the 
u(N f )n c Knizhnik-Zamolodchikov equations, which in the replica limit (Np — > 0) are given 
by equations (|A22|) with k = Nc- For example, setting all replica indices equal to 1 yields 

(g li (i)g t 1I (2)g| i (3)g li (4)) ~ |T| 2 [k Nc (z)k Nc (i -z) + k Nc (i - z)k Nc {z)\ (3.10) 

where T = [zuZ23 z (1 — z)Y l ^ N c and Kn c {z) and En c {z) are natural generalizations of 
the complete elliptic integrals — see equation (|A19| ). In addition one may consider 'mixed' 
correlation functions of the form: 

(g 1I (i)g t li (2)gt 5 (3)g 25 (4)) ~ |T| 2 [e No { z )(k Nc {i -z)- e No (i - z)) + c.c.) (3.11) 

where c.c stands for complex conjugation — replacement of z by z. Correlation functions cal- 
culated with respect to this effective (replicated) action correspond to correlation functions 
averaged with respect to the initial action with quenched disorder p!")[|Ti| : 

(Q(l)Qt(2)Qt(3)Q(4)) ^ = (g rir - i (l)g^(2)g^- i (3)Q rir - i (4)) rep (3.12a) 
(Q(l)Qt(2)) A (Qt( 3 )Q(4)) A = (Q nr - (l)Qt ifi (2)Qt. r -.(3)Q rjrl (4)) rep , n ± Tj (3.12b) 
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The left hand side of ( |3.12a|) represents the four-point correlation function of the Q-field 
(calculated with respect to the original action in the presence of quenched disorder) averaged 
over disorder (denoted by an overline); the left hand side of ( |3.12b| ) represents the product 
of two two-point correlation function of the Q-field (again calculated with respect to the 
original action in the presence of quenched disorder) averaged over disorder. 

We note that as a direct consequence of the rather simple replica index structure of ( p.9|) 
the results ( |3.12a| ) and ( |3.12b| ) are independent of which replicas are actually considered; 
this is a manifestation of replica symmetry — see for example §3.3 of reference [p5 ]. The 
correlation functions of the LDOS may be obtained from these results by means of the 
decomposition ( |3.1| ) together with crossing symmetry. We emphasize that one does not 
perform a trace over replica indices in order to extract the LDOS; the traces over replica 
indices appearing in equation (60) of reference |31| and equation (5) of reference [^3| are 
erroneous. We further draw attention to the simplicity and manifest crossing symmetry of 
the results ( |3.10| ) and ( |3.11| ) as compared with those obtained in reference 



B. Strong Disorder Approach 



As was first discussed in p4|,p9|| , and subsequently developed in [pf| , the random Dirac 
fermion problem is amenable to a direct treatment in the limit of infinite disorder strength 
(gA — > oo) without invoking either replicas or supersymmetry. We note that in this limit 
the probability measure ( |2.3| ) is absent. We summarize here only the important details. We 
separate the action fl2TT|) into its chiral components 



, N C 

S= d 2 iY, l R U B + iA ) R <* + L «(d + iA)La\ 

J a=l 



(3.13) 



and parameterize the gauge fields in terms of group elements g(£) G SU C (N) ~ SL(N; C) 
residing in the complex extension of SU(N) 



A = idgg \ A = idgg 1 . 
This enables one to perform the chiral gauge transformations 

L-> gC, gTZ, 



(3.14) 



(3.15) 



so as to render a theory of free fermions (7Z, C) decoupled from the gauge fields. An im- 
portant feature of this procedure is that the Jacobian of the transformations ( |3.15| ) is pro- 
portional to the partition function of the original action ( |2.1| ) at fixed disorder, Z[A^\ 
this cancels the normalizing partition function in (fixed disorder) correlation functions and 
removes the need to invoke replicas or supersymmetry in order to perform disorder averag- 
ing. The Jacobian associated with the change of variables ( |3.14j ) is well known to involve 
the WZNW action (§!]) on the (non-compact) manifold h = g^g e SU C (N)/SU(N) at level 
k = 



2 A c 34,39,561, 



VA = Vg exp (2N c W- 2No tf g]) 



(3.16) 
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In reference |37j the Wakimoto free- field representation of the SU C (N)/ SU(N) WZNW 
model at k — —2N was constructed, and it was shown to generate the su(N)__ 2 n Kac- 
Moody algebra; in other words, the correlation functions of the /i-fields may be obtained 
from the solution of the su(N)_ 2 n Knizhnik-Zamolodchikov equations — see appendix |A 2 



In this approach the Q-fields of the random Dirac fermion problem with Nc colours are 
expressed as primary fields of the su(Nc) -2N C WZNW model 'dressed' by free fermions: 

N c N c 
a,a=l a,a=l 

We shall study the conformal dimensions and correlation functions of these fields in the 
subsections below. 



1. Conformal Dimensions 



It is readily seen from the strong disorder decomposition (|3.17| ) that the conformal di- 



mension of the Q-field is that of free Dirac fermion {h = 1/2) and an su(Nc)-2N c primary 
field — see equation (|A4| ) and set N = Nc and k = —2Nq: 



1 N c - 1 1 



Hq ~ 2 + 2N C (N C - 2N C ) ~ 2N C (3 " 18) 



This agrees with the replica result ( |3.8|) and ensures the equality of the two-point and 
three-point correlation functions in both approaches. The coincidence of higher correlation 
functions will be discussed below. 



2. Correlation functions of the Q-field 

The correlation functions of the Q-field are thus obtained by a fermionic 'dressing' of 
the correlation functions of the su{Nc)-2N c WZNW model. We note that whilst the chiral 



solutions to the su(Nc)-2N c WZNW model given in [37]] are correct, their normalization 
constants are erroneous — they do not satisfy the conformal bootstrap. We provide the 
correct conformal blocks in equations ( |A25|) and their expression in terms of generalized 



elliptic integrals in equations (|A26|). Moreover, in the light of the new results obtained 



in the replica approach, we shall find it convenient to generalize the decomposition (|3.17|) 



slightly so as to accommodate disorder averages of products of quenched correlation functions 
such as those appearing in equation (|3.12b| ). To this end we introduce as many additional 
fermionic species (denoted by an index p) as quenched correlation functions we wish to 
disorder average. We emphasize that these additional indices are not required to perform 
the disorder averaging (as would be true of replicas) but simply encode which combinations 
of quenched correlation functions are to be averaged over disorder — in this approach we 
do not perform an N p — > limit. The Q-matrix acquires a pseudo-replica index structure 
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(labelled by the index p) analogous to that appearing in equation (|3[ 

N C N C 

Qp,p ^ ^Q:,P ^a,& Qp,p ^ ] £ , a,ph'a,a'^'tt,pi (3.19) 

a,a=l a,a=l 

The correlation functions of the fields (|3.19|) may be written in the form 

2 

(Q PlP - 1 (l)Q P2P - 2 (2)QL 3 (3)Q P4 p 4 (4)) = \z u z 23 r 4h £ f#>(*) (3.20) 

a,b=l 

where h — (1 — N C )/2N C is the conformal dimension of the su(Nc)-2N c primary field h a ^ 
as given by equation ( |A4| ) . Gp a \zi) is the (holomorphic) product of the free Dirac and 
su(Nc)-2N c correlation functions traced over the colour indices: 

=£ ( — - — ) X (l?F a \z) + I?F a \z)) (3.21) 



a=l 



The coefficients C a f, appearing in equation (|3.20|) have the values Cn = C21 = 1, C\\ = C22 = 
to ensure single- valuedness and their off-diagonal form reflects the logarithmic nature of the 
underlying su(N c )-2N c LCFT. We use the symbol p here to denote the ordered sequence 
of pseudo-replica indices Pi,P2,P3,Pi and the invariant tensors I\ and I 2 are defined as 
If = S Pl ' P2 6 P3,p ' 1 , 1 2 = S Pl,P3 S P2 ' P4 ; analogous expressions hold for a. Performing the trace 
over colour indices appearing in (|3.21|) by means of the identities 

Nc N C 

J^I?I? = N 2 C J^I?I« = N C (3.22) 

a=l a=l 

and utilizing the relations ( |A29|) satisfied by the su(Nc)-2N c chiral blocks one obtains: 

Q?{zi) = ^—A{l P E Nc (z) + I%[K Nc (z) - E Nc (z)] } (3.23a) 

^14^32 

= — a{j? \e Nc {z) - K Nc (z)} - I p 2 E Nc (z)\ (3.23b) 

^14^32 L J J 

where d = cN(N 2 — 1) is a normalization constant, A = [z(l — z)]~ l l N c , and where we have 
adopted the notation that f(z) = f(l — z) for an arbitrary function f(z). In particular 
by inserting the explicit results ( |3.23j ) into the decomposition ( |3.20| ) and collecting the 
coefficients of I P I P one may recast (|3.20f) so as to read: 

2r' 2 2 

where the Fij(z, z) are (both fortunately and remarkably) single- valued combinations of the 
solutions to the u(0)n c Knizhnik-Zamolodchikov equations as given by equations 



with k = Nc- That is to say (upto an irrelevant normalization) we have recovered the 
replica result ( |3.9| ) in which our psuedo-replica indices play the role of replica indices. In 
particular the results ( |3.10| ) and (|3.11|) together with their interpretations ( |3.12|) follow 
straightforwardly. 
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C. Supersymmetric Approach 



In the supersymmetric approach to disordered fermionic systems one introduces bosonic 



copies of the original Grassmann fields [42, 43]. The partition function of the resulting 



supersymmetric theory is equal to unity due to the inverse relationship between ordinary 
c-number Gaussian functional integrals and their Grassmann counterparts. The absence of a 
disorder dependent partition function normalizing quenched correlation functions drastically 
simplifies their disorder averaging. The supersymmetric approach to the random Dirac 
fermion problem has been outlined by Bernard and LeClair J58| who, following the general 



principles of the supersymmetric approach, have introduced bosonic copies of the Grassmann 
fields coupled to the same gauge potential. The resulting action is given by (|3.3|) where the 
summation over replicas is replaced by a summation over fermionic (Grassmann) and bosonic 
(c-number) fields: 

N c 

gSUSY 



d 2 iY, E r'P a V (3-25) 

„-_1 „, /3_1 



i=l a,/3=l 



The symmetry of the free supersymmetrized Dirac action (in the absence of any gauge 
fields) is OSp{2N c \2N c ) H and it may be recast as the 6sp(2N c \2N c )i WZNW model. In 
particular, the random su(N c ) gauge potential couples only to currents J a (and J a ) residing 
in the su(Nc)o subalgebra of the complete 6sp{2Nc\2Nc)\ Kac-Moody algebra: 

5 SUSY = osp{2N c \2N c )i + J d 2 ^ (J a A a + J a A a ) (3.26) 

In the special case Nc = 2, Bernard and LeClair have demonstrated that the Suguwara 
energy momentum tensor for osp(2Nc\2Nc)\ may be decomposed into the sum of two com- 
muting pieces pertaining to different symmetries 



^~osp(4|4)i — 2~osp(2|2)_ 2 + ^su(2) (3.27) 

Exploiting the decomposition (|3.27|) one may rewrite equation (|3.26|) in the decoupled form 



S SVSY = 55p(2|2)_ 2 + su{2) + / d 2 i {J a A a + J a A a ) (3.28) 

We note that the role of the decomposition ( |3.27| ) in this approach closely mirrors that of 
the decomposition ( |3.6|) employed in the replica approach — both allow us to decouple the 
effects of the gauge potential disorder. Indeed, the latter two terms appearing in equation 
( [3.28j ) are precisely those which appear in the replica approach as Np — > 0. As was rigorously 
proven in references []3U],|3T|, and utilized to our advantage in section [111 A| , this su{2) sector 



becomes massive for the simple Gaussian distribution of the gauge fields given in equation 
( [2.3p ; a perturbative renormalization group argument was given in reference [^] where 
the one loop beta function was calculated with the result that the coupling constant 
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flows to strong coupling.^] That is to say, the low-energy effective theory governing the Q- 
field is the osp(2|2)_ 2 WZNW model. Fortunately, the 6sp(2\2) k WZNW model has been 
discussed quite extensively in the work of Maassarani and Serban p9| |. In particular the 
model undergoes a dramatic simplification at k — —2 |K| and we provide a rather extensive 



discussion of this model in Appendix |B[ In this approach the Q-field may be represented as 

Q = Q 1 ' 1 , g 1 " = g 4 ' 3 (3.29) 



where the Q a ' 
see Appendix 



are primary fields transforming in the [0, 1/2] representaion of osp(2\2) 



1. Conformal Dimensions 

In the supersymmetric approach to the random su{2) Dirac fermion problem, the con- 
formal dimension of the g-field coincides with that of a primary field transforming in the 
fundamental [0, 1/2] representation of <5ip(2|2)_ 2 — see equation flET| ) and set k = —2: 



h Q = - = - (3.30) 

y 4-2-2 8 v ; 



This is in agreement with the replica result (|3.8|) and the strong disorder result ( p. 18 ) when 
N c = 2. 



2. Correlation functions of the Q-field 

The four-point correlation function of the g-matrix admits the osp(2\2) x osp{2\2) in- 
variant decomposition — see equations (fB4l), ( [Bo]) and flB8l): 

3 

(Q a ^{l)Q a ^{2)Q a ^{3)Q a ^{A)) = \z 14 z 23 \- 1/2 ^ IflfF^z) (3.31) 

where a denotes the ordered sequence of indices a 1 ,a 2 ,« 3 ,a 4 , which label the basis states 
of the four-dimensional representation of osp(2\2), and where the invariant tensors I\, J 2 
and ^3 are defined in equations (p9| ) together with similar equations for Ji, J 2 and I3. The 
anharmonic ratio z is defined as z = z^z^/ z 14 z 32 and similarly for z. The functions Fij(z, z) 
are single- valued combinations of the solutions to the o?p(2|2)_ 2 Knizhnik-Zamolodchikov 



Although a formal proof of the gap generation exists only for the case when the disorder distri- 
bution is Gaussian, it appears reasonable to assume that the su(2) sector remains massive for a 
broader choice of disorder distributions. At scales smaller than the gap, the effective action is given 
by the critical osp(2|2)_ 2 WZNW model. The critical point is stable with respect to variations of 
the disorder: all such variations affect only the massive modes and hence do not generate relevant 
perturbations of the critical action. 
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equations, which are summarized by equations (|B29| ) — (|B33|) . For example, focusing on 
the correlation function pertaining to the decomposition Q3.29Q : 

(Q 1I (l)g 43 (2)Q 43 (3)Q 1I (4)) ~ |T| 2 [K{z)K{\ - z) + K(l - z)K(z)] . (3.32) 

where here T = [2142232(1 — z)] -1 / 4 and K(z) and E(z) are the complete elliptic integrals. 
This is in agreement with both the replica result ( p.lOj ) and the strong disorder result ( 3.24 ) 
when the number of colours Nc = 2. 



IV. DENSE POLYMERS AND TWIST OPERATORS 



It turns out that the supersymmetric description provides an extremely economical and 
straightforward approach to the disordered Dirac fermion problem. It follows from the work 
of Rasmussen fl6l| , that the action of the <53p(2|2)_ 2 WZNW model may be represented as 



the direct sum of three simple theories: 



o^(2|2)_ 2 = S«(2) 1 + i- J d 2 i(d ll i P ) 2 + J d 2 £e ab d, X a d»X b 



(4.1) 



where <p is a non-compact bosonic field, x a is a two-component symplectic fermion fl62 |, 
and e a fe is a two-component antisymmetric tensor. As expected in this supersymmetric 
theory, the bosonic sector (c = 2) and the symplectic fermions (c = —2) together yield 
a total central charge of zero. We note that the presence of the su(2)i WZNW model in 
the decomposition ( f4.1|) of the <55p(2|2)_ 2 WZNW model reflects an underlying su(2)-k/2 
Kac-Moody subalgebra residing in the 6sp(2\2)k algebra [^|. The decomposition ( f4.1|) may 



be simplified even further by noting that the su(2)i WZNW model admits the following free 
field representation: 



^(2) 1 = i- f d 2 ad. 



(4.2) 



where is a compact free boson; note that our choice of normalization is for latter conve- 
nience in equation (|4.9| ). That is to say, the oSp(2|2)_ 2 WZNW model may be represented 
as the sum of a compact bosonic field (c = 1) a non-compact bosonic field (p (c = 1) and 
a two-component symplectic fermion (c = —2): 



55p(2|2)_ 2 = / d 2 £ 



7^0) 2 + ^ d ^f + ta b d, X a d»X b 



(4.3) 



In view of the decomposition C N=.3|) one anticipates a representation of the osp(2|2)_ 2 primary 
fields in terms of the primary fields of the models appearing on the right-hand side of equation 
( |4.3|) — namely vertex operators e a ^ and and the primary fields of the c = —2 theory. 
As we shall discover in sections IV A - IV (J this is indeed possible. The non-unitary c = —2 
minimal model has received a great deal of attention in recent years as a theory of dense 
polymers [p3|,B5 , as a celebrated example of a logarithmic conformal field theory [|6^ ,|66|-|70 



and as a conformal ghost system ||71|| . The structure of this theory is rather rich and is 
known to consist of several sectors. As we shall see below, the so-called Z 2 twisted sector 
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of the c = —2 theory will play a crucial role here; this sector consists of a scalar primary 
field fi of conformal dimension —1/8 and tensor primary fields of conformal dimension 
3/8. The correlation functions and operator product expansions of these fields have been 
studied in reference p4| and more fully in reference |69|]. In section |IV A| we shall recast 
the four-point correlation functions of the o3p(2|2)_2 WZNW in a form which facilitates 
the elucidation of the desired operator correspondence. In section |IV B| we shall discuss the 
Z2 twist operator correlation functions of the c = —2 model [5D|. In section |1V we shall 
compare the correlation functions (summarized in Table | and Table |I| of sections [IV A| and 
IV B| respectively) and arrive at the aforementioned correspondence. 



A. The oip(2|2)_ 2 WZNW model 

In view of the large number of components of the generic <5sp(2|2)_ 2 non-chiral four-point 
correlation function and the rather cumbersome and opaque invariant tensors ( |B9| ) we shall 
refashion the correlator somewhat.0 In particular the c = —2 operator correspondence will 
follow quite naturally. We study the four-point correlation function of the supersymmetric 
Q-matrix: 

F a '*( Zi ,Zi) = (g Ql ' Sl (^l^l)Q Q2 ' (i2 (^,^)Q Ql ' (i3 ^3^3)Q a4 ' Ci4 (^^4)), (4.4) 

where on the left hand side we use the symbol a to denote the ordered sequence of indices 
«i, a 2 , a 3 , a 4 . The indices assume the values 1,2,3,4 and label the basis states of the 
four- dimensional [0, 1/2] representation of osp(2\2); in the notation of [BUJ states 1 and 4 are 
bosonic, whilst states 2 and 3 are fermionic. As may be seen from our more detailed studies 
in Appendix [BJ, this correlation function may be written in the (off-diagonal) form 

J^ a (zi,Zi) = -^ 1 \z l )^ {2 \z{)-^ 2 \z t )^ 1 \z l ), (4.5) 

in which we have set the overall normalization of the four-point function to minus unity for 
latter convenience, and where 

3 

= Iz^r^^/fif 1 ^). (4.6) 
i=i 

Using the explicit form of the chiral blocks F^ a \z) appearing in equations (B2"3), ( B2(j| ) and 
( |B27| ) together with the known transformation properties of the invariant tensors — see 
equations flB48|) and ( |B50|) — we deduce that 

F*>( 2 \z i ) = -VF k '< ? -\z i ). (4.7) 



2 Restricting our attention to four-point correlation functions pertaining to the four-dimensional 
[0, 1/2] representation we have a total of 4 8 components. Since the vast majority of these vanish 
however, and many are related by symmetry, it is desirable to distill them further. 
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Sector 


a 


v v 


JDOSOIllC 


1 1 A A 

1144 


A A 1 1 

441 1 


yid^f ) 1 ZJ\ \Z ) 


1414 


4141 


(Apry,-! T (1 — A K(r\ 


1441 


4114 


-(4e 7 )- 1 Tir(^) 


Fermionic 


2233 


3322 


(4e7j 1 \ltj\z) - {Z- z)K (z)\ 


2323 


3232 


(4e 7 ) T [2E(z) - (1 - z)K(z)} 


2332 


3223 


(4e7j l \irj\z) - A [z)\ 


Mixed 


1234 


1324 


2143 


3142 


±T[E{z)-{l-z)K{z)\ 


2413 


3412 


4231 


4321 


1243 


1342 


2134 


3124 


±T[K{z)-E{z)} 


2431 


3421 


4213 


4312 


1423 


1432 


2314 


3214 


±TE{z) 


2341 


3241 


4123 


4132 



TABLE I. Chiral correlation functions of the 53p(2|2)_2 WZNW model displaying regu- 
lar behaviour in the vicinity of z = and logarithmic behaviour in the vicinity of z = 1; 
T = [Z14Z23 z(l — z)] -1 / 4 and z = 212-234/2142:32 — see Appendix |B| for further details. 



We use the tilde to denote the interchange of the coordinates or indices 2 and 3, and V 
to denote the fermionic parity of this permutation; note that this permutation induces the 
transformation z — > 1 — z in these functions. Substituting equation ( |4.7| ) into equation ( |4.5| ) 
we obtain the result 

T a '«{z h Zi) =VF a > { - 1 \z i )F~ & ^ 1 \t i ) +VT^ 1 \z i )^' { - 1 \z i ). (4.8) 

That is to say, one may build the full quota of single-valued and crossing symmetric non- 
chiral correlation functions of the osp(2\2)_ 2 WZNW from our knowledge of the chiral func- 
tions ^"'(^(zj) which display regular behaviour in the vicinity of z = and logarithmic 
behaviour in the vicinity of z — 1. Using the explict results for the chiral blocks (|B23|), 



( P26| ) and ( P27| ) together with the invariant tensors (|B9|) we are able to gather these non- 



trivial functions in Table |. We emphasize that the functions J 7 "'^^) are not simply the 
conformal blocks of the o?p(2|2)_2 WZNW model, but also embody the non-trivial tensorial 
structure of osp{2\2) — see equation 



B. Twist Operator Correlation Functions 

The correlation functions of the twist operators /i and v in the non-unitary c = —2 
minimal model have been extensively studied by Gaberdiel and Kausch [|55 |. Their non- 



chiral four-point functions are built from linear combinations of their chiral counterparts 
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Sector 


Correlator 


Regular contribution 


n 








\ IS IS IS IS I 


T 3 z [2_Ef,zl — C2 — zliffVll 


V 




T 3 (1 - *) [2E(z) - (1 - z)lf(z)] 




\ IS IS IS IS 1 


T 3 [2ftYzl - K(z)\ 




(/iZ/ ± Z/ T /i) X Z M 1/2 Z 2 + 3 1/2 
/,/±„,„/¥\ v -1/2 


±T[E{z)-{l-z)K{z)\ 


Mixed 


{^^) x Zn /2 z£ /2 
(z^/izA/z) x z 13 z 24 


±T - 




( / u/i^ ± ^ T ) x z^ l2 z^ 12 

(U^fifi) X z£ /2 z£ /2 


±T£(z) 



TABLE II. Chiral correlation functions of the Z 2 twist operators of c = —2 displaying reg- 
ular behaviour in the vicinity of z = and logarithmic behaviour in the vicinity of z = 1; 
T = [Z14Z23 z{\ — z)\~ l / A and z = Z12Z34/Z14Z32 — see Gaberdiel and Kausch for further details [69] 
and note that our definition of the anharmonic ratio differs from theirs. 



so as to respect the stringent constraints of single-valuedness and crossing symmetry. As 
may be seen from §4 of their work, the contributions from either chiral sector are typically 
composed of two elliptic integral solutions posessing logarithmic branch cuts; one of these 
solutions displays regular behaviour in the vicinity of z = and logarithmic behaviour in 
the vicinity of z — 1, whilst the other displays regular behaviour in the vicinity of z — 1 
and logarithmic behaviour in the vicinity of z = 0. As in section [IV Aj > in order to establish 
our correspondence Q4.9| ) we find it convenient to focus on the chiral contributions to the 
four-point functions which exhibit regular behaviour in the vicinity of z — 0. The chiral 
four-point functions most relevant to our discussion and displaying such regular behaviour 
in the vicinity of the origin are summarized in table 



C. Twist Operator Correspondence 

Comparing our correlation functions for the osp(2|2)_2 WZNW model (summarized in 
Table [I]) with those of Gaberdiel and Kausch for the Z2 twist fields of the c = — 2 minimal 
model |69[ (summarized in Table |TJ) we see a clear correspondence at the level of their el- 
liptic integrals. In particular we see that the bosonic and fermionic states of the dsp(2|2)_ 2 
WZNW model are naturally associated with the Z 2 twist fields // and v a respectively. The 
remaining powers of Zy which distinguish Table | from Table |Tj are provided by bosonic ver- 
tex operators as suggested by the decomposition (|4.3|). Explicitly, we establish the following 
equivalence between the chiral operators of the oip(2|2)_ 2 WZNW model and the 'dressed' 
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chiral operators of the c = 



—2 minimal model:^ 



Q 1 ~A~ 1 e i +ti, Q 2 ~Ae-' p iy-, Q 3 ~Ae v v + , Q 4 ~ A' 1 n, (4.9) 

where A = (4e7) 1//4 is a free parameter corresponding to the arbitrary relative normaliza- 
tions of the su(2) doublet (Qi, Q4) and the two singlets (Q2, Q3) in the four- dimensional 
representation of osp(2\2) [pH|, [72|. With the normalization chosen in equation ( |4.3| ) the 
compact bosonic exponents e ^ have conformal dimension = 1/4, and the non-compact 
bosonic exponents e 93 have /i^ = —1/4. It is thus straightforward to see that the dimensions 
add up to the correct value of 1/8 — equation ([B7|) — in the decomposition ( f4.9|) . 



We have now highlighted the important role played by the c = —2 model in the chiral 
structure of the osp(2|2)_ 2 WZNW model, and have provided a prescription for construct- 
ing the non-chiral correlation functions (|4.8|) . Although there are many interesting and 
notable exceptions, we emphasize that in general, the non-chiral correlation functions of the 
oSp(2|2)_2 WZNW are not simply related to single non-chiral correlation functions of the 
Z2 twist operators dressed by non-chiral bosons. For example, it follows from our solution 
of the osp(2\2)_2 Knizhnik-Zamolodchikov equations that 

(Q 2 ' 2 (1)Q 3 ' S (2)Q 2 - 2 (3)Q 3 ' S (4)> = C{[(1 - z)K - 2E]{{1 + z)K - 2E] + c.c.}, (4.10) 



where C = (16e 2 7 2 ) 2 1 214 223] ^ 2 A, and A is given by equation (|B33|) . We note that this is 



not simply related to a single correlation function of the non-chiral operators fj,(z, z) and 
v a ,a( z , 2). That is to say, no single non-chiral four-point function of the Z 2 twist operators 
gives rise to this particular single-valued combination of the elliptic integrals. We invite the 



reader to verify this statement with the aid of §4 of Gaberdiel and Kausch [B9|. It may be 



interesting to study the non-chiral aspects of this correspondence in greater detail. 



V. CONVERGENCE OF APPROACHES 

In section |l| we demonstrated that three different approaches to the random Dirac 
fermion problem yield identical results for the four-point correlation functions of the local 



density of states. In section |V| we have further demonstrated that the supersymmetric ap- 
proach inherits its non-trivial logarithmic structure from the c = —2 non-unitary minimal 
model. In this section we shall discuss the convergence of the these approaches — summa- 
rized in table [D I - and the emergence of the c = —2 theory in a little more detail. Indeed, 
we shall point to a remarkable proliferation of c = —2 theories and effects. In table |ITl| we 
have highlighted the separation of the active (critical WZNW) degrees of freedom from the 
passive (massive or decoupled Jacobian) degrees of freedom in each approach. 

An interesting aspect of this marriage of approaches is the agreement between the weak 
coupling (replica and supersymmetry) limit and the strong coupling limit. Although such an 



3 We use the (overused and rather abused) symbol ~ to emphasize that we have derived the 
correspondence from the chiral regular solutions and that it is valid up to phase. 
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Approach 


Fermions 


Bosons 


Active 


Passive 


Replicas 


N c x N F 





u(l) + su(N F ) Nc 


su{N c )n f 


Strong Disorder 


N c 





+ su{N c )i + su(N c )-2N c 


(c = -2)^-1 


Supersymmetry 


2 


2 


6?p(2 1 2) _ 2 


sm(2) 



TABLE III. Current algebra approaches used to investigate the critical behaviour of the random 
non-Abelian Dirac fermion problem and the separation of the active (critical WZNW) degrees of 
freedom from the passive (massive or decoupled Jacobian) degrees of freedom. The central charges 
and the conformal dimensions may be seen to add up correctly in each approach. 



agreement might have been anticipated purely on physical grounds, its emergence is quite 
remarkable from a field-theoretic perspective. At the very least, it is quite surprising that 
such superficially different conformal field theories may nonetheless yield equivalent results, 
albeit for a specific subset of physically motivated correlation functions. In particular, the 
resultant central charge of the active degrees of freedom in the strong disorder approach 
differs from that in the replica and supersymmetric approaches, namely zero. As we have 
discussed in section [Tl [B| — and indicate in table |I! I - the active degrees of freedom in the 



strong disorder limit may be expressed as the sum of three different models: 

S = u{l)+su{Nc)i + su{Nc)- 2 N c - (5.1) 

where the first two models describe Nc colours of free massless Dirac fermions, and the 
remaining model encodes the non-trivial logarithmic structure. We use the well known 



result for the central charge of the g^ WZNW model |53 

k dimg 



(5.2) 



where g v is the dual Coxeter number of the algebra g (equal to N for su(N)) together with 
the fact that a free boson u(l) has central charge 1, to find the resultant central charge of 

1 + N c -2N C + N C ° ° r 1 ; 

That is to say, the active degrees of freedom in the strong disorder approach yield a positive 
central charge, in contrast to the replica and supersymmetric degrees of freedom which yield 
a net zero central charge. 

As we have discused in section [TV], the supersymmetric approach to the disordered Dirac 
fermion problem reveals a hidden substructure which is inherited from the c = —2 minimal 
model. In view of the convergence of approaches outlined in this paper it follows quite 
naturally that the active degrees of freedom in the replica and strong disorder treatments 
also inherit non-trivial traits from the c = —2 model. In addition to the rather natural 
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proliferation of elliptic integrals in the conformal blocks presented in this paper, we note 
that a relation between the sw(2)_ 4 WZNW model (which arises in the strong disorder 
treatment with Nc = 2) and the c = —2 theory has already been discussed in the string 



theory literature | 73[| : the sw(2)_ 4 WZNW model is cohomologically equivalent to an N = 4 
supersymmetric bosonic string with c = —2 matter.^ 

A curious by-product of these investigations is that the presence of the c = —2 theory 
may apparantly be seen in the 'passive' degrees of freedom outlined in table |TTT[ Indeed, a 
number of c = —2 models (arising as Jacobians) explicitly decouple in the strong disorder 
approach and play a passive (or spectator) role in regards the LDOS. A possible explanation 
for the closely related structure of the active and passive theories indicated in table |TJ may 
arise from the requirement of the mutual cancellation of their logarithmic singularities |60| . 
In particular we anticipate a relationship between the sw(2)o WZNW model and the c = —2 
theory and a correspondence analogous to equation ( |4.9|) . We note that the four-point 
functions of the su(2) model have been studied in references [T75|,[77[] and assume a simple 
form in terms of the complete elliptic integrals The investigation of the su(2) WZNW 
model will be continued in reference [7S| . 

As we close this section we comment briefly on the stability of the critical point with 
respect to variations of the disorder distribution. In the replica and supersymmetric ap- 
proaches the gauge potential disorder is coupled directly to the massive sector and it is natu- 
ral to assume that the critical theory is protected from such variations. In the strong disorder 
approach, however, the disorder variations are coupled directly to the critical su(Nc)~2N c 
subsector, and its stability is far from obvious. It is therefore interesting to study the sta- 
bility properties of this critical theory. A simple perturbation of the critical theory that one 
might consider is the deformation of the WZNW action by its kinetic term. This operator 
has scaling dimension zero and is strongly relevant. However, at level k = —2N (or more 
generally —2g v ) this operator commutes with all the Kac-Moody currents and therefore 
does not affect the correlation functions — see Appendix 6 of reference |T2| and note the 
different sign conventions for the level k. In a more general framework, the WZNW models 
at level k = —2g v arise quite naturally in both two-dimensional (c = 0) topological field 
theories possessing a non-Abelian current algebra |79j , and in the strong disorder treatment 
of arbitrary WZNW models coupled to random vector potentials [39|| . The critical level 
arises from a BRST (Becchi-Rouet-Stora-Tyutin) |§0|,§TJ symmetry nilpotency condition 
and is required for the coexistence of a Kac-Moody algebra symmetry and a topological 
algebra symmetry J7^] . We emphasize that the level of the underlying Kac-Moody algebra, 
the existence of a topological algebra and the stability of our strongly disordered theory are 
therefore intimately related. 



4 We note that this matter + string theory admits a description in terms of the c = — 2 model 
dressed by c = 28 Liouville theory and c = — 26 string ghosts. We remind the interested reader 
that this (and other) Liouville LCFTs f|,|74|,|75| emerged in a very natural way in the closely 
related problem of prelocalization in disordered conductors [35,35]. 
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VI. CONCLUSIONS 



We briefly summarize here a number of our main results: 

• We have solved the random non-Abelian Dirac fermion problem by means of three 
different non-perturbative procedures based on the replica approach, supersymmetry, 
and in the limit of strong disorder. We have demonstrated that this menage a trois of 
approaches yields identical results (for the four-point correlation functions of the local 
density of states) in this relatively simple, but quite non-trivial model of disorder. 

• We have emphasized the special role played by the level k of the Kac-Moody algebra. 
As we have seen both here and in reference [KJ the level k = —2 is rather special for the 
osp(2\2)k WZNW model which arises in the supersymmetric treatment of the random 
Dirac fermion problem: an entire conformal block decouples from the spectrum and 
the resulting theory is drastically simplified. 

• We have established that the c = —2 minimal model plays an important role in the 
random non-Abelian Dirac fermion problem. We have found a rather simple and 
suggestive form for the supersymmetric critical action ( f4.3|) : 



oSp(2|2)_ 2 = it(l) + gl{\) + [c = -2]. 

We have highlighted the relevance of the twist operators of the c = — 2 theory in 
the disorder averaged correlation functions of the Dirac fermion problem. The c = —2 
symplectic fermions are ubiquitous in the construction of supersymmetric sigma models 
[ft2l |43|l; this example indicates that this sector may well be responsible for many 
remarkable features of disordered critical points including the presence of logarithmic 
operators ||66|| . We draw attention to the fact that the c = —2 theory emerges in the 



theoretical description of dense polymers [63],E5|. We hope that this connection may 



ultimately yield a more intuitive picture of the random non-Abelian Dirac fermion 
problem and other disordered critical points. 

We have argued for the stability of the critical point with respect to variations in 
the disorder, and noted that the level of the underlying Kac-Moody algebra, the 
existence of a topological algebra and the stability of our strongly disordered theory are 
rather intimately linked. It is an interesting open question whether such considerations 
may yield Kac-Moody level selection mechanisms in models displaying lines of critical 
points [0 . 
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Note added in proof: whilst this paper was being completed we became aware of a 



preprint by A. W. W. Ludwig |82j in which one of our results, namely the equivalence ([O] 
was established. 
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APPENDIX A: THE «u(JV) fc WZNW MODEL 



In this appendix we consider the appearance of logarithms in the correlation functions of 
the su(N)k WZNW model. Following §4 of Knizhnik and Zamolodchikov []53| we compute 
the four-point functions 



J 7a ' a (zi, Zi) = (g(z 1 ,z 1 )g\z 2 ,z 2 )g\z 3 ,z 3 )g(z 4 ,z 4 )). i 



(Al) 



of the field g(zi,Zi) = g ai,ai (zi, Zi) transforming in the fundamental representation of 
SU(N) x SU(N). We use the symbol a. to denote the ordered sequence of indices aii, a 2 , a 3 , 014. 
Global conformal invariance restricts this correlation function to have the form 



(z 14 z 23 z 14 z 23 )- 2h F a ^(z,z) 



where z and z are the anharmonic ratios 

_ £12 £34 . 

Z — , A 

^14^32 

and the conformal dimension h of the field g is 

A 2 



^12^34 

^14^32 



h 



The correlation function 



2A(A + k) 

admits the su(A) xsu(A) invariant decomposition 

2 

^(z,z) = Y,I?I?Fv(^z) 

ij=l 



(A2) 



(A3) 



(A4) 



where the invariant tensors I\ and I 2 are defined as 



(A5) 



(A6) 



with similar equations for I\ and I 2 . The four scalar functions satisfy the coupled first- 



order differential equations [E3 



dF_ 

dz 



z z — 1 



Q 



where the matrices P and Q are given by 

1 



P 



N(N + k) 



N 2 - 1 N 
-1 



F, where F 



Q 



Vj 







N{N + k) \N N z -1 



(A7) 



(A8) 



There are similar equations for the antiholomorphic dependence. Suppressing the antiholo- 
morphic index j from the functions F±j and F 2 j, one may obtain the second-order differential 
equation satisfied by F\(z): 

N 2 (N + k) 2 z 2 (l- z) 2 F['(z) - 

N(N + k)z(l -z)[2- N(2N + k) + (3A 2 + Nk — 4)z\ F[{z) + (A9) 
[1 - N 2 - (A 4 - 6 A 2 — Nk + A)z + (A 4 - 5 A 2 + A)z 2 ] F x {z) = 0. 
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The corresponding F 2 (z) m &y be obtained from the F±(z) solutions: 



F 2 (z) = -(N + k)zF[(z) + N~ x (l - z)- 1 [l - N 2 + (N 2 - 2)z] F^z). (A10) 

Upon the change of variables 

F 1 {z) = z- N ^ N+k \z{l - z^^G^z) (All) 

one obtains an equation of the hypergeometric form: 

z(l - z)G'[ + [7 - (a + j3 + l)z] G[ - apGi = 0. (A12) 

where a = —1/(N + k), (3 = 1/(N + k) and 7 = k/(N + k). For non-integer values of 7 
the solutions of these equations are given by Knizhnik and Zamolodchikov — see equations 
(4.10a) and (4.10b) of p3| . For integer values of 7, however, the solution involves logarithms 



the roots of the indicial equation corresponding to ( A12|) (namely and 1 —7) differ by an 



integer. Of particular importance in the study of the disordered Dirac fermion problem are 
the cases — * (7 = 1) appearing in the replica treatment (section [111 A| ) and k = —2N 



(7 = 2) appearing in the strong disorder treatment (section |IIIB|) . We consider both of 
these cases below. 



1. The u(N) k WZNW Model as N 



As may be seen from equation (|A4j) the conformal dimensions of the 5u(A^)^ WZNW 
model diverge as A^ — > 0, as do some of the prefactors in the chiral blocks ( |A11| ); the A^ — » 



limit of the su(N)k WZNW model alone fails to yield a well defined CFT. However, the 
replica approach to the random Dirac fermion problem instructs us to consider the u(N)k 
WZNW model. As we shall demonstrate below, the u(N)k WZNW model does have a well 
defined A^ — > limit. 

An arbitrary element u, of the group U(N), may be expressed as an element of SU(N), g, 
multiplied by a phase: u = e taip g. Using the Polyakov-Wiegmann identity for the WZNW 
model |R3| 



W k [ab] = W k [a) + W k [b) + — I d 2 £Ti'(a- 1 dabdb- 1 ) (A13) 

2ix J 

with a = e lOLLp and b = g one obtains^ 

n 2 Nk f 

W k [u] = W k [g] + —— / d 2 i d^ip (A14) 



5 We have used the facts that the current gdg 1 residing in the su(N) algebra is traceless, thereby 
eliminating the second term in ( |A13| ), that the WZNW term ( p.^| ) vanishes for the U(l) element 
e iaip , and that Tr' = 2Tr. 
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The u(N)k WZNW model is therefore seen to be the su(N)k WZNW model augmented by 
a free scalar field of prescribed normalization. The conformal dimension of the field e taip is 
fixedby this normalization to beQ/i a = l/(2Nk). The conformal dimension of the composite 
field u = e iaip g is thus 

1 N 2 -1 1 

hu(N h = ^ + m{N + k) ^ ^ (A15) 

and is seen to have the finite replica limit l/2k 2 . In addition, the (non- vanishing) holomor- 
phic four-point correlation function of the U(l) fields V a = e iaip is given by 

(V a (l)V_ tt (2)V- a (3)Va(4)) = {z liZ , 2 )-^ Nk [z{l - z)]~ 1/Nk (A16) 

The divergences occuring in the SU(N)fc results flA2j ) and (|A1 1|) as N — > are seen to 
be compensated by those of the U(l) phase ( |A16| ). In the replica (N — > 0) limit the 
hypergeometric equation ( |A12|) reads 

z(l - z)G'{ + (1 - z)G[ + k- 2 d = (A17) 

and one is able to find the N — > limit of the u(N)k conformal blocks: 

if } = 7 2Fl [-H;l;z} (A18a) 
F[ 2) = 1 {l-z) 2 iq[l-± l + i;2;l-z] (Al8b) 



fc> 1 fe' 

^ (1) =7| 2 Fi[l-il + i;2;z] (A18c) 



Ff = 7 fc 2 F 1 [-i,i;l;l-z] (A18d) 
where 7 = [z(l — z)]^ 1 ^ 2 . We find it convenient to introduce generalizations of the complete 



elliptic integrals of the first and second kind [84]: 



^jEp^l-ii;^] (A19a) 
E k (z) = l 2 F 1 [-l,l;l;z] (A19b) 



In terms of these functions (and rescaling) the N — > u(N)k conformal blocks read: 

^ = 7^(2) (A20a) 

if > = 7 fc [A^(l - z) - £ fc (l - z)] (A20b) 
7 [AT fc (z) (A20c) 

F {2) = 1 kE k {l-z) (A20d) 



6 For a free scalar field governed by the action S = \g J d 2 ^ d^ipd^cp, the field e iaip has h a = a 2 /8irg 
— see for example equations (5.73) and (6.60) of [l7fl. 
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The full U (N) x U (N) invariant correlation function is built from single- valued combinations 
of these conformal blocks; it is straightforward to see that one can only have 



F^z, z) = X 12 F^(z)F 2 \z) + Fl 2 \z)F 1 \z) 



Substituting the explicit forms (|A20|) into ( |A21|) one obtains 



F ll = X 12 k\ 1 \ 2 
F l2 = X 12 k\ 1 \ 2 
F 2 i = X 12 k \>y\ 2 
F 22 = X 12 k\ 7 \ 2 



E k (K k -E k ) + E k {K k -E k ] 
E k E k + {K k -E k ){K k -E k ) 
E k E k + (K k -E k )(K k -E k ] 
E k {K k -E k ) + E k {K k -E k ) 



(A21) 

(A22a) 
(A22b) 
(A22c) 
(A22d) 



where we have adopted the notation that f(z) = f{l — z), f(z) = f(z) and f(z) = f(l — z) 
for an arbitrary function f. 



2. The su(N)- 2N WZNW Model 



As may be seen from equations ( |A1 1| ) and ( |A12| ) with k = —2N } the Fi (z) chiral blocks 
for the su(N)^ 2N WZNW Model are given by 

F l {z) = z[z{l-z)]- l/N2 G 1 (z) (A23) 
where G\ satisfies the hypergeometric equation with a = 1/N, (3 = —1/N, 7 = 2: 

z{l - z)G'[ + (2 - z)G\ + N~ 2 G 1 = (A24) 



Solving this equation and applying the relation ( |A10|) enables one to obtain the full set of 
su(N)_ 2N chiral blocks:[] 



Az 2 Fi[-j; i;2;z] 



F[ 2) = A (1 - z) 2 2 iq[l - 1 1 + 1 ; 3; 1 - 2 



,(2) 



Az 2 
'~2N 



Fl [l-j,,l + j,;3;z 



-2NA(l-z) ^ 1 i;2;l-« 



(A25a) 
(A25b) 

(A25c) 

(A25d) 



7 We note that ( A/25b ) may be obtained by substituting G\{z) = z~ l (l — z) 2 H\{z) into (|A24 ) and 
replacing z by 1 — z — this yields a hypergeometric equation. The result ( ]A25cf ) may be obtained 
by straightforward application of ( A10| ) to ( A25a| ). In deriving ( |A25d ) one may utilize the result 
(N 2 -l)z(l-z) 2 iq[2-i 2 + ^;4;z]+3iV 2 (2-z) 2 F 1 [1 - ^ 1 + ^; 3; z] = 6iV 2 2 F X [ 
which is easily verified by means of the Gauss recursion relations. 



_ J_- 

" TV' N 



;2;z 
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where A = [z(l — z)\ . In terms of generalized elliptic integrals the su(N)_ 2 n chiral 
blocks read: 



Fi 1] = cA [[1 + (N — l)z] E N -(1- z)K_ 
2NcA \[1 + (N — l)z] E N - NzK_ 



(2) 



,(1) 
2 

,(2) 



-cA [[N -(N- l)z] E N - N(l - z)K N ] 



-2NcA 



[N — (N — l)z] E N - zk 



N 



(A26a) 
(A26b) 

(A26c) 
(A26d) 



where c = N 2 /\n[N 2 — 1)]. We have used the Gauss recursion formula GR 9.137(13) p5| 

(A27) 



7[a - (7 - f3)z\ 2 F x [a,f3;y, z] - 07(1 - z) 2 F x \a + l,/3;7; z] 
+ (7 -a) (j-/3)z 2 F 1 [a,(3; 1 + l;z] =0 



with a = -1/N, (3 = 1/N, 7 = 1 to obtain (|A26aD , and GR 9.137(6) gj 

7(7 + 1) 2 F, [a, P; r,z] -7(7 + 1) 2 F t [a, 13] 7 + I; z] 
-a(3z 2 F 1 [a+ 1,(3 + 1] 7 + 2; z) =0 



(A28) 



with a = -1/N, (3 = 1/N, 7 = 1 to obtain (^26^ ). One may obtain flA~26TJ) and ( |A"26c[D 
by noting that F 2 (1) (^) = -F x (2) (l - z)/2N and F 2 {2 \z) = -2A^F 1 (1) (1 - z) as follows from 
( |A25p . In particular we note that 

N 2 F[ l) + iVF 2 (1) 
NF[ l) + N 2 F 2 (1) 
N 2 F^ + NF 2 (2) 



c'A^^jv 

c'k{l-z)[K N -E N \ 
-2Nc'Az[K N - E N ] 
-2Nc'A(l - z)E N 



(A29a) 
(A29b) 
(A29c) 
(A29d) 



where d = cN(N 2 — 1) is simply another constant. 



APPENDIX B: THE osp(2|2) fc WZNW MODEL 

1. Representations of osp(2|2) 

The Lie superalgebra osp(2|2) ~ spl(2|l) consists of four bosonic generators Q3, Q+, Q-, 
B, and four fermionic generators W+, W-, V + , V— The bosonic subalgebra is si (2) © u(l): 

[Q 3 , Q±] = ±Q±, [Q + , Q_] = 2Q 3 , [5, Q±] = 0, [5, Q 3 ] = (Bl) 

and the eigenvalues of Q3 and S (called isospin and baryon number respectively) are used 
to classify the basis states of finite-dimensional representations [f72], 186]. A representation 
[b, q] contains at most four multiplets of states: 

\b,q), |6+i,g-|), |6-|,g-i), \b,q-l), (B2) 
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the states within a given multiplet \b, q) being labelled by their third component of isospin 
\b,q,q 3 ). In particular the [0, |] representation is four-dimensional and contains a doublet 
and two singlets (the multiplet \b,q — l) being absent): 

[0, |, |>, [0, |, — 1>, ||,0,0), |-|,0,0). (B3) 



Following ]55[] we denote these states as |1), |4), |3) and |2) respectively. In addition to their 
isospin and baryon number, the basis states carry a grading denoted by e = ±1. States |1) 
and |4) are even/bosonic and have e = 0, whilst states |2) and |3) are odd/fermionic and 
have e — 1. On this basis of states one may construct 4x4 matrix representation [0, ^] of 
the generators of osp(2|2). 



2. Knizhhnik Zamolodchikov Equations 

Let us study the four-point function of the supersymmetric Q-matrix: 

F^(z u Zi) = (g Ql ' Sl (^l^l)Q a2 '' i2 (^^2)Q ai ' S3 (^3^3)Q a4 ' <54 (^4^4)), (B4) 

where on the left hand side we use the symbol a to denote the ordered sequence of indices 
ati, «2, «3, «4- Global conformal invariance restricts this correlation function to have the 
form 

F^izi, Zl ) = (z u z 23 z u z 23 )- 2h F a ^(z } z) (B5) 
where z and z are the anharmonic ratios 

Z12Z34 _ ^12^34 / T~\ \ 

z = , z = _ (Bo) 

^14^32 ^14^32 

and the conformal dimension h of the field <fr is 

h = T^Tk < B7 » 

The correlation function (|B5|) has the osp(2|2) xosp(2|2) invariant decomposition 

3 

F<*>*(z,z) = Y,I?I?FiAz,z) (B8) 

ij=l 



The tensors I and I are given in appendix A of ||59|| : 
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h = (1144) + (1234)4e7 + (1324)4e7 - (1414) + (2143)4e7 

+ (2233)16e 2 7 2 + (2323)16eV - (2413)4e7 + (3142)4e7 + (3232)16eV 

+ (3322)16e 2 7 2 - (3412)4e 7 - (4141) - (4231)4e 7 - (4321)4e 7 + (4411) (B9a) 



I 2 = (1234)4e7 - (1243)4e7 + (1324)4e7 - (1342)4e7 

- (1414) + (1441) - (2134)4e 7 + (2143)4e 7 + (2233)32e 2 7 2 + (2323)16e 2 7 2 

+ (2332)16e 2 7 2 - (2413)4e7 + (2431)4e7 - (3124)4e7 + (3142)4e7 

+ (3223)16e 2 7 2 + (3232)16eV + (3322)32eV - (3412)4e7 + (3421)4e7 

+ (4114) - (4141) + (4213)4e 7 - (4231)4e7 + (4312)4e7 - (4321)4e 7 (B9b) 



h = (1234) - (1243) + (1324) - (1342) + (1423) + (1432) - (2134) + (2143) 
+ (2233)8e 7 + (2314) + (2323)8e7 + (2332)8e 7 - (2341) - (2413) + (2431) 

- (3124) + (3142) + (3214) + (3223)8e7 + (3232)8e7 - (3241) + (3322)8e 7 

- (3412) + (3421) - (4123) - (4132) + (4213) - (4231) + (4312) - (4321) 

and the nine scalar functions i*y satisfy the coupled first-order differential equations 



(B9c) 



x 



dF 

dz 



P 



z-1 



Q 



where F 



7 3„ 



Vj 



(BIO) 



There are similar equations for the antiholomorphic dependence. Suppressing the anti- 
holomorphic index j from the functions F±j ■ ■ ■ F 3 j , one may reduce this first-order matrix 



differential equation to the following set of equations [59 



x 3 z 3 (l 



•d 3 F 3 (z) 



+ x 2 (l + 2x)z 2 (l - z) 2 (l - 2z)^fi- + 



dz 3 

xz(l — z)[— 1 



X 



2xz 



dz 2 

. , .,dFs(z) 
2x(2 + x)z(l - z)]—^ + 



F*(z) 



1 



Ae^xz{l — z) 



[-l-x + 2z + 2xz{l - z)]F 3 (z 
[x 2 D 2 F 3 (z) + 2x(l - z)DF 3 (z) + (1 - 2z)F 3 (z 



Fi{z) = — [xDF 3 (z) - F 3 (z)} + (z- 2)F 2 (z) 



(Blla) 

(Bllb) 
(Bile) 



where D = z(l — z)d/dz. Equation (pila|) has three independent solutions which we denote 



(a) 



F 3 where a = 1,2,3. Equations (pilbj ) and (P3 1 1 c| ) yield the corresponding solutions F, 
and F x (a) . The nine scalar functions Fij appearing in equation (|B8|) may be expressed as a 
linear combination of these nine functions (the so-called chiral/current blocks.) 



F lJ (z,z)= Y,X ab Ft\z)Ff\z 



(B12) 



a,b=l 



The values of the coefficients X a b are determined by single valuedness (monodromy invari- 
ance) and crossing symmetry to be discussed later. 
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a. Factorization at Level k = 1 



At level k = 1 (x — 1) the Knizhnik-Zamolodchikov equation ( [Blla] ) takes the form 

z 3 (l - z)F% + 3z 2 (l - 2z)Fg - 2z(l + 3z)Fz - 2F 3 = (B13) 
This may be factorized to read 

P (2) D (1) F 3 = (B14) 

where 



X?« = ^ + 1 

P( 2 ) =z 2 (l-z 



d 2 



d 



(B15) 

3z2 dz- 2 (B16) 
It is straightforward to see that the equation D^'F^ = 0, and thus (|B14j ) has the solution 



(B17) 



Applying the formulae ( |Bllb| ) and ( |Bllc[ ) to this single solution of the Knizhnik- 
Zamolodchikov equation (P14 ) one obtains: 



2z 



Aejz(l — z) 



(i) 



z-2 



4ej(l - z) 



(B18) 



Analysis of the full set of solutions of ( |B14| ) reveals that it is only these blocks which enter the 
physical correlator - a fact which is intimately related to the factorization of the Knizhhnik- 
Zamolodchikov equation (|B14 ). In the su(2) WZNW model the factorization properties of 



the differential operators are related to a finite closure of the underlying operator algebra [87 



b. Factorization at Level k 



At level k = —2 (x = 4) the differential equation ( |Blla| ) takes the form 

[4z(l - z)fF'z + 9[4z(l - z)} 2 (l - 2z)Fl + 

4z(l - z )[-5 + 8z- 48z(l -z)]F' 3 + (B19) 

[-5 + 2z + 8z(l - z))F 3 = 

This equation may be factorized to read 

D«p(2)_F 3 = o (B20) 

where this time the differential operators T>^ and T>^ take the form 
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£> (1) = 4z(l - z)^ - (5 - 6z) (B21) 

P( 2 ) = [4z(l - z)] 2 ^ + 8z(l -z){3- Az)^- + (1 - 4z) (B22) 

In the light of the factorization which occured at k = 1 one might expect that the subset 
of solutions obtained from V^F^ = 0, which are closed under analytic continuation on the 
Riemann sphere, might serve as a reduced basis on which to perform the conformal bootstrap. 
This indeed turns out to be the case. It is easily seen that the equation V^F 3 = 0, and 
therefore the full Knizhnik-Zamolodchikov equation ( |B2U| ) admits the two solutions^] 

F m (A _ E(z) (2) _ E(l-z)-K(l-z) 

F * {z) -[ Z (i-z w 3 () ~ — w^w 1 — ( } 

where K{z) is the complete elliptic integral of the first kind and E(z) is the complete elliptic 
integral of the second kind|] 



1 ^ 



zx 



2 



E(z)= / dx K(z)= / dx (B24) 

Jo Jo y/(l-X 2 )(l-ZX 2 ) 

The representation (|B23| ) is particularly useful owing to the very simple manner in which 
the elliptic integrals behave under differentiation with respect to the parameter z: 



dE(z) E{z) - K{z) dK{z) E{z) - (1 - z)K(z) 



dz 2z dz 2z(l — z) 

Applying ( |Bllb| ) and ( |Bllc| ) to these solutions yields: 



(B25) 



,(i) _ K(z) (i) _ zK(z) 



2 ~" 4e 7 [, 2 (i-, 2 )] 1 /4 1 '" 4e 7 [z(l-z)] 1 /4 ^ 
(2) = K(l-z) (2) = zK(l - z) 

2 4e 7 [z(l - z)} 1 / 4 1 4e 7 [z(l-z)] 1 /4 1 

As we shall subsequently demonstrate, one may satisfy the demands of single-valuedness 
and crossing symmetry on the subspace of functions ( |B23| ), ( |B26| ) and ( [B27| ). Once again 
this is intimately connected with the factorization of the Knizhnik-Zamolodchikov equation 
(B20). Having found closed form expressions for the chiral blocks, one must now construct 



8 Upon the change of variables F3 = [z(l — z)] 1 / 4 H(z), the equation V^F^ = reduces to 
the canonical form of the hypergeometric equation z(l — z)H + [c — (a + b + l)z]H — abH = 
with a = -1/2, b = 1/2, c = I. This has solutions = 2 iq [-\,\;l\z] and F (2) = (1 - 
z) 2F1 [5, |; 2; 1 — z] . Using the well known fact that 2-^1 [— \, 2, 1; z] = ^E(z) together with it's 
derivative, the result follows. 

9 Note that many texts on the theory of elliptic integrals denote the parameter z by k 2 - the 
so-called modulus. This is purely a mater of convention. 
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the physical correlation functions to be single-valued on the whole Riemann sphere. It is 
enough to ensure this property at the two singular points z = and z = 1. The blocks (z) 
are regular at z = and logarithmic at z = 1, whilst the blocks F> (z) are logarithmic at 
z = and regular at z = 1. It is straightforward to see that in this subspace of functions 
one can only have[^ 



F ij (z,z) = X 



12 



(B28) 



One may gather the explicit expressions for the Fij at level k = —2 in the Hermitian matrix: 



\z\ 2 F 2 2 -zF 22 -zF- d 



23 



-zF 22 
-zF 23 



F 22 

F 2 3 



F 2 3 

F33 



(B29) 



in which we have singled out the elements 



F 22 = -A 
F 23 = Ae-fA 



KK + KK 
KR + KE — KE 



F 33 = (4e 7 ) 2 A 



EE - EK + EE - EK 



(B30) 
(B31) 
(B32) 



where 



A 



X 



12 



(4e 7 ) 2 |^(l -z)|V2 



(B33) 



and where we have adopted the notation that f(z) = /(l — 2), /(z) = /(z) and /(z) 
/': I — z) for an arbitrary function /. One may also demonstrate that this combination is 
consistent with the crossing symmetry constraints on the correlator: 



Z,Z) 



WF & > & {1 -z,l- z), F a ' a {z, z) = z~ 2A z- 2A VVF & ^(l/z,l/z) (B34) 



where a. denotes the permuted sequence of indices a 1; 03, a 2 , a±, et denotes the sequence 
ai, a 4 , a 3 , a 2 , and V denotes the parity of the permutation: 



V 



V 



(B35) 



The proof of this requires the use of the analytic continuation formulae of the elliptic integrals 
and is presented in appendix |B~1. 



10 A detailed proof of this statement in terms of monodromy matrices is given in appendix B 3 
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3. Monodromy Invar iance 

A monodromy transformation of a function of z consists in letting z circulate around 
some other point (typically a singular point). We define 



C F(z,z) = lim F(ze 2iw \ze- 2i7Tt ) (B36) 
t-»i- 

& F(z, z) = Urn F(l + (z- \)e 2m \ 1 + (z - l)e- 2iwt ) (B37) 

Using the standard analytic continuation formulae for the hypergeometric series, it is easily 
seen that the elliptic integrals have the following nontrivial monodromy properties: 

C K(1 -z) = K(l -z)- 2iK(z) (B38) 

C E(1 -z) = E(l -z) + 2i[E(z) - K(z)} (B39) 

&K(z) = K(z) - 2iK{\ - z) (B40) 

dE(z) = E(z) + 2i[E(l -z)- K{\ - z)\ (B41) 

together with the trivial transformations CqK(z) = K(z), CqE(z) = E(z), C\K(1 — z) = 
K{1 — z), CiE{l — z) = E(l — z). Using these results it is straightforward to see that 

C F l {a \z) = (g ) ah F t (b \z), i = 1,2,3 (B42) 
C 1 Ft\z) = (g 1 ) ab F?\z), i = l,2,3 (B43) 

where (on this reduced subspace) the matrices go and g\ are given by 

»-(1ty< »-(o < B «> 

Under the monodromy transformation Co, the combination 

Fy(z,z)= j^X ah Fl a \z)Ff\z) (B45) 

a, 6=1 

transforms in the following manner 

CoF^z, z) = F^(z)F^(z) [I n - 2t(X 12 - X 21 ) + 4X 22 ] + F^\z)Ff\z) [X, 



22 



F?\z)F? \z) [X 12 + 2iX 22 ] + F?\z)F?\z) [X 21 - 2iX, 
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Invariance under the monodromy transformation Co thus requires X\ 2 = X 2 ±, and X 22 = 0. 
That is to say 

F lJ (z,z)=X 11 F! 1 \z)Fj 1 \-z)+X 12 [fP(z)fP(z) + F^(z)F^(z) 
Under the monodromy transformation C\ this simplified function transforms as 

KFyiz, z) = F?\z)F?\z) [In] + F^\z)f\ 2 \z) [X 12 - 2iX u ] 
+ F?\z)fU(z) [X 12 + 2tX n ] + F?\z)Ff\z) [4X n ] 



31 



Invariance under the monodromy transformation C\ therefore imposes the additional con- 



straint X\\ = 0. Hence monodromy invariance restricts F t 



F ij (z,z)=X 12 F^(z)F^(z)+Fr(z)F^'(z) 



(2), 



>J 



<2). 



.Z,Z. 



to have the form 



(B46) 



as stated in the text. 



4. Crossing Symmetry 

a. Invariance under z — > 1 — z 

Crossing symmetry requires that 

F a > & (z, z) = VVF^ h {\ - z ,l-z) (B47) 

where a denotes the sequence of indices ai,a 2 ,a 3 , a 4 , a. denotes the permuted sequence of 
indices ai, a 3 , a 2 , a 4 , V = (— l) £a 2 £a 3 denotes the parity of the interchange in the holomorphic 
sector (e a is for bosons and 1 for fermions.) Introducing the following tensor |59|] 

J? = Vlt (B48) 
the crossing symmetry constraint may be written 

3 3 

ItlfFijiz, z) = J?J?W - M - *) (B49) 

i,j=l i,j=l 

The tensor J admits the following decomposition 

J« = C\ 3 If d = f -1 1 -4e 7 ] (B50) 

Substituting this decomposition into equation ( |B49| ) and equating the coefficients of I^Ij on 

both sides, one finds the following nine identities which must be satisfied by the Fij(z,z) if 

this crossing symmetry is to be satisfied. Denoting Fij(z, z) by and Fij(l — z, 1 — z) by 
FiA these follows: 

F u = F n + F 12 + F 21 + F 22 (B51a) 

F 12 = -F 12 - F 22 (B51b) 

F 13 = 4e 7 (F 12 + F 22 ) + F 13 + F 23 (B51c) 

F 21 = -F 21 - F 22 (B51d) 

F 22 = F 22 (B51e) 

F 23 = -Ae-iF 22 - F 23 (B51f) 

F 31 = 4e 7 (F 21 + F 22 ) + F 3 i + F 32 (B51g) 

F 32 = -Ae-fF 22 - F 32 (B51h) 

F 33 = 4e7(4e 7 F 22 + F 23 + F 32 ) + F 33 (B51i) 

It is straightforward to show that these relations are indeed satisfied. 
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b. Invariance under z 
Crossing symmetry requires that 

>> z) = 



-2A; 



(B52) 



where a denotes the sequence of indices aj, «2, «3, «4, ck denotes the permuted sequence 
of indices ai,a^, «3,«2, P = (— l) eQ 2( e a3+ £a 4)+ e a3 e a 4 denotes the parity of the interchange in 
the holomorphic sector (e a is for bosons and 1 for fermions.) Introducing the following 
tensor 

(B53) 



the crossing symmetry constraint may be written 



1/X--1/X 



(B54) 



The tensor admits the following decomposition |59|j 



Co 




(B55) 



Substituting this decomposition into equation (P54|) and equating the coefficients of Iilj on 



both sides, one finds the following nine identities which must be satisfied by the F^iz, z) if 
this crossing symmetry is to be satisfied. Denoting Fij(z, z) by F^, and Fij(l/z, 1/z) by F^ 
these follows: 

(B56a) 
(B56b) 
(B56c) 
(B56d) 
(B56e) 
(B56f) 
(B56g) 
(B56h) 



F u = 


\z 


- 2/x F 2 2 






F\2 = 


\z 


-*/»F n 






F 13 = 


\z 


- 2 ^[-4e 7 (F 21 


+ F22) 


— F23] 


F21 = 


\z 


~ 2/x F 12 






F22 = 


\z 


-2 /x p n 






F23 = 


\z 


- 2 ^[-4e 7 (Ai 


+ F 12 ) 


-A3] 


F31 = 


\z 


- 2 ^[-4e7(A 2 


+ F22) 


— F32} 


F32 = 


\z 


-2/*[_ 4e7 (^ 1 


+ F 21 ) 


-F 31 ] 


F33 = 


\z 


- 2/ "[16eV(Ai + F 12 


+ F21 + 



4e7(F 13 + F 23 + F 3 i + F32) + F< 



33J 



(B56i) 



In order to demonstrate that these identities are satisfied we shall make use of the following 
rather simple transformation laws of the elliptic integrals under z — > 1/ 



We note that replacing i by —i in equations ( |B57a ) and (B57b) changes the domain of validity 
of the transformation from Jmz < to 3mz > 0. Since E and K appear together in the Fij it is 
essential that their transformations be defined in the same region of the complex plane. 
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K(l/z) 
E{l/z) 
K(l - l/z) 
E(l - l/z) 



z +1/2 [K{z) + %K{\ - z)} , 3mz < 
[D(z) - iD(l - z)} , 3mz < 



z -l/2 

Z +V*K(1 - z) 
z-^EiX-z) 



where 



D(z) = E{z) - (l-z)K{z). 



(B57a) 
(B57b) 
(B57c) 
(B57d) 



(B58) 



These are easily obtained using the standard analytic continuation formulae for ordinary 
hypergeometric functions. Recalling the form of F 22 appearing in equation (B30| ) namely 



F 22 = -A [K(l - z)K{z) + K{z)K{l - z)\ 
it is straightforward to see how it transforms under z l/z: 



22 



22 



(B59) 



(B60) 



Since Fu = \z\ 2 F 22 from equation ( |B29| ) it follows that that constraint ( [B56a| ) is satisfied. 
Replacing now z by l/z in ( [B56a|) one may infer the validity of ( |B56e| ). Further, since 
F 2 i = —zF 22 from equation ( |B29|) it follows by reciprocity that 



21 



-zr x F 22 . 



Substituting for F 22 using ( |B60[ ) and using the fact that F 12 = —zF 22 one obtains 



21 



1 1/2 p 



12 



(B61) 



(B62) 



That is to say, constraint ( |B56b| ) is satisfied. By reciprocity we see that constraint (|B56d| ) 
is also satisfied. Recalling now the form of F 2 % appearing in equation (P31|) namely 



F 23 = Ae-fA KK + KE - KE 
it is straightforward to see how it transforms under z — > l/z: 
p 23 = 4e T k| 3/2 A [z 1 / 2 {k + iK) z 1/2 R+ 

z^kz- 1 ' 2 (d 



(B63) 



iD 



l l 2 K + ik)z- l l 2 E 



4e 7 |2| 1/2 A 

4e 7 |^| 1/2 A 
4e 7 |^| 1/2 A 



K + ik)K + zk(E-(l- z)K) 



+izK [E -zK) - z[K + iK) E 



I 2 ( KK + KK) +z(kE- KK - KE) 



KK + KK 



z (KK + KK )+z[ KK + KE-KE 



4e 7 |^ 



1/2 



-\z\ 2 F 22 + zF 22 + —F 23 



z 
4e 7 " 



I 1 / 2 [-4e 7 (Fn + F 12 ) - F, 



13 



(B64) 

(B65) 
(B66) 
(B67) 

(B68) 
(B69) 
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Replacing z by I/2 on both sides of this equation one obtains the relation (|B561j) . Taking 
now the complex conjugate of ( |B56f|) and using the Hermiticity property = Fj iy one 
obtains ( |B56h| ). Further, rearranging ( |B69| ) for F\ 3 and replacing F\\ and F\ 2 using ( |B56a ) 



and ( |B56b| ) respectively one obtains (|B56c|) . Taking now the complex conjugate of ( |B56c 
and using Hermiticity, one obtains (|B56g|). We consider finally how F33, namely 



(B70) 



F 33 = (4e 7 ) 2 A E(E -K) + c.c. 

behaves under the transformation z — > 1/z: 

F 33 = (4e 7 ) 2 |z| 3 / 2 A \z-V 2 (D - ib){z-^ 2 E - z^K) + c.c. 
= (4e7) 2 |z| 1/2 A \{D - iD)(E - zK) + c.c. 



(4e7) 2 |z 
(4e 7 ) 2 |z 



V2 A 
V2 A 



EE — (1 — z)KE - zEK + z(l - z)KK + c.c. 

1 - z\' z KK + c.c. 



EE - EK + (1 - z){EK -KE + KK) 

l ' 2 [(4e 7 ) 2 |l - z\ 2 F 22 + 4e 7 [(1 - z)F 23 + c.c] + F 33 ] 

l ' 2 [(4e 7 ) 2 (F n + F12 + F 2 \ + F 22 ) + 4e 7 (F 13 + F 23 + F 31 + F 32 



+ F 33 ] 



(B71) 
(B72) 
(B73) 
(B74) 

(B75) 
(B76) 



which may be seen to be the relation (|B56i| ) with z replaced by 1/z. This completes our 
proof of crossing symmetry in the reduced subspace. 
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